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ABSTRACT 

We suggest a technique for construction of profiles of incoherent pulsar emission using 
instantaneous angular distribution indicatrix of radiation from relativistic source moving 
in the neutron star magnetosphere. In the general case the indicatrix depends on the 
kinematic parameters of radiation source such as velocity and acceleration. The method 
is illustrated by calculation of the profiles of pulsar radiation based on synchrotron 
radiation indicatrix. This technique can be easily generalized for the case of other types 
of relativistic particles radiation. The considered profiles are compared with observed 
ones. A good agreement was found with profiles of 7-pulsars — Geminga, Crab and 
Vela. 



1. Introduction 

There is no clear conception of the nature of pulsar electromagnetic radiation. We accept the 
well known point of view suggested by I.S. Shklovsky (Shklovsky) that at least the high frequency 
part of radiation represents the incoherent synchrotron radiation. 

There are different models of pulsar radiation based on various physical processes which are 
sometimes contradictory and insufficiently conditioned. We suggest a technique of kinematic pro- 
jecting of synchrotron radiation (Bagrov 1965) which gives the pulsar radiation profiles includ- 
ing the phase dependence of polarization characteristics (Bagrov et al. 1972; Bagrov et al. 1974; 
Bordovitsyn). Then one can analyze the coherence effects and possible mechanisms of pulsar radi- 
ation. 

The suggested technique of determination of the pulsar radiation profiles is as a matter of fact 
a kinematic projecting of relativistic radiation regardless to the specific source of the radiation. The 
role of such sources of radiation can play relativistic charged particles electrons, magnetic particles 
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(neutrons) or even bunches of those particles. The main source of information in this technique 
is the highly directed indicatrix of the angular distribution of instantaneous power of relativistic 
radiation, which is specified by the kinematic and dynamic characteristics of the emitting object 
in every specific case. It is known (Shklovsky; Manchester and Taylor) that there are two models 
of radiation geometry. The first one is due to two beams of curvature radiation associated with 
the magnetic poles. The second one is associated with synchrotron radiation in a plane orthogonal 
to the magnetic axis of pulsar. We restrict our consideration with the second model of radiation, 
though it is possible to apply this technique in the first case too. 



2. Transformation of angular distribution into the pulse profile 

We start from the indicatrix of synchrotron radiation averaged over the period. Because of high 
energy main part of radiation is concentrated near the plane of symmetry which we call the light 
plane. Directions, assigned to radiation maximum, are situated in the light plane. Axis orthogonal 
to the light plane coincides with the direction of pulsar magnetic moment fj,. Vector fx is precessing 
around the pulsar axis of rotation M with the angular velocity Q. 

We assume that the angular velocity u> of the emitting particles is much more than CI. 

eH 

L0= >fi, (1) 



here e is the charge, mo is the mass of rest, c is the speed of light, 7 = yl — f3 2 is the relativistic 
Lorentz factor, (5 = u/c, u is the particle velocity. It follows from the last inequality that 

eH 

7 « o- (2) 

Let us estimate the order of 7-factor satisfying the last condition. It can be taken for pulsars 
that Q ps 10 s _1 , H ps 10 12 Oe not far from the surface and H ss 10 6 Oe in the periphery of 
magnetosphere (near the light cylinder). Then the condition 2 gives the value 7 <C 10 12 -7- 10 18 , i.e. 
condition 2 is always satisfied in available relativistic and even ultrarelativistic case. 

We introduce the coordinate system shown in Fig. 1. Direction of the angular velocity of pulsar 
is parallel to Z-axis. The pulsar magnetic axis defined by unit vector k is specified by the angles k 
and fit. The unit vector of the observer's direction is given by the angles 5 and A. We assume that 
vector k is situated in the plane YZ at the moment t = 0. 

We present the angular distribution of radiation averaged over the period of particle motion 
in a form 

dW = W p{9) sm6d9, 

where W is the total power of radiation, which is relativistic invariant, p{9) is the function of the 
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angular distribution of radiation normalized to unity 

7T 

J p{6) shxdde = 1. 



In the case of synchrotron radiation (Bagrov et al. 1972) 

W = ~ — w 2 /?V = W SR 



and 



P{0) 



o 7 



2 + (3 2 sin 2 9 

;i - (3 cos ef/ 2 



sin 2 9 [4 + (3 2 sin 2 0) 1 
4 7 2 (1 -/3 cos 0) 7 /2 J 



(3) 



It is obvious that the angle theta varies in the course of time, as the unit vector k setting the 
position of light plane is precessing around the Z axis. It is easy to show that 



cos9(t) = cos 5 cos k + sin S sin k s'm(£lt — A). 



(4) 



If we know 9(t), we determine the radiation profile pit) and the power dW emitted in the direc- 
tions of the angle interval sin 9 d9. The maximum of radiation will be observed at the moment of 
intersection of the observation direction by the light plane when cos 9 = 0, i.e. when 



sin(J7t — A) = — cot 5 cot k. 
The latter equation has a solution when | cot 5 cot k\ < 1. 



(5) 



Equation (5) has two roots within the precession period which are equal at 5 = ir/2 ± k. In 
other words the light plane intersects the observation direction twice within the period, if the angle 
between the vector n and plane XY is less than k. 

If the light plane intersects the observation direction (or touches it), then the angle A can be 
chosen such that the observation direction has been situated in the light plane at t = 0. Under this 




Fig. 1. — The coordinate system. 
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initial condition the radiation pulse takes on the maximum value p(9 = tt/2) at t = 0. Thereat it 
follows from the equation (5) that 

sin A = cot 5 cot re. (6) 

Another form of relation (6) is possible 

cos 5 , 
sin/t = , < re, o < tt/2. 

V cos 2 5 + sin 2 5 sin 2 A 

It follows from equation (6)that for the same acute angles re and S < tt/2 — re the term sin A 
takes equal values within the range < A < tt, i.e. Ai and A2 = tt — X\, which corresponds to two 
possible intersections of the observation direction and the light plane. At Ai = A2 = tt/2 we have 
re = tt/2 — 5. In this case the observation direction only touches the light plane. Note that equation 
(6) is symmetric relative to angles re and S. It means that the positions of the two peaks of profiles 
will not change if we swap re and 5. The case of re = is special. It follows from equation (4) that 
in this case cos = cos 5 =const and pulsar is observed as a source of synchrotron radiation with 
constant luminosity and angular distribution of radiation defined by relation (3). 

In the ultrarelativistic case of 7 S> 1 the main part of radiation is concentrated near the light 
plane 9 ~ tt/2. Thereat duration of the radiation pulse is much less than f2 _1 . It allows to expand 
the function p(8) in a series in small angles ip = tt/2 — 0, or small Qt. To a first approximation in 
we get 

ijj = cos 5 cos re + sin S sin re(f2i cos A — sin A) . (7) 
Specifying the initial conditions in accordance with equation (6), we have 

tp = Qt sin 5 sin re cos A, (8) 

or, with account of (6), 

ip = ±Qt\/ sin 2 5 sin 2 re — cos 2 6 cos 2 re. 
We introduce another angular variable 

^ = jtp = ~{ut, uj = f2sin<5 sin re sin A. (9) 

The variable tJ has a sense of angular velocity at which the observer's direction intersects the light 
plane. The highest velocity corresponds to passage through the radiation maximum. In this case S 
and A are associated with the direction in the light plane. 

Substituting equations (8) and (9) into equation (3) we find the radiation profile in a form 

3 7+12 X 2 

P{t) = V 7 (l + x 2 )^ - (10) 

If(5 = 7r/2±reor^ = re then uJ = and we have to take into account the terms of higher order in 
the expansion in fit. 
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Figs 2 display the synchrotron radiation profiles correspond to weakly relativistic electron at 
(3 = 0.9 and at various values of angle 5 and fixed k. Fig. 2 a shows the case when the angle k = 90° 
i.e. when vector k is orthogonal to rotation axis. If the observer is situated in the direction of Z 
axis (<5 = 90°), radiation of constant intensity is registered. As the angle 5 increases the radiation is 
modulated at the frequency of 2VL and the effect of glimmer pulsar is observed. One can see that the 
greater is the angle 5 the deeper is the glimmer modulation. If the angle between pulsar magnetic 
axis and rotation axis is acute (see Fig. 2 b), then the distance between the peaks of radiation within 
a period decreases with increasing of angle S, and at 5 + k = tt/2 they merge in one maximum. At 
further decreasing of 5 the light plane does not intersect the observer's direction and the amplitude 
of splashes falls quickly off. When the angle between the pulsar magnetic axis and rotation axis 
is small, the light plane intersects the observation direction under a small angle (see Fig. 2 c). In 
consequence the splashes duration increases considerably and the depth of modulation falls off. 
Thereby the main part of radiation is concentrated near the plane orthogonal to the rotation axis 
(<5 « 90°). 

Fig. 3 shows the radiation profiles of ultrarelativistic particles {[3 = 0.99). They are distin- 
guished by deeper modulation and small pulse duration. In Fig. 4 the light profiles calculated by 
use of approximated equation (10) (solid line) and exact equation (3) with account of (4) (dashed 
line) are shown for particles moving at the velocity (3 = 0.999. We see that even for moderately 
relativistic particle these equation are in a good agreement. 

3. Comparison with the observed profiles 

We have considered a simple model of pulsar emission. Some of the obtained profiles are close 
in shape to really observed profiles. In this section we compare the profiles of 7-ray pulsars, the 
Crab, Vela and Geminga with the profiles obtained from synchrotron radiation indicatrix. As one 
can see from Figs. 3 and 4 the variation of angles k and S gives a wide variety of profiles. But it 
does not mean that you can produce any profile in this way. In particular, profiles obtained by 
this method have a certain relationship between the width of peaks, phase position and the level of 
plato between them. Thus, not every observed profile can be built with help of used model. But if 
a calculated profile fits an observed one, it can mean that the suggested model is adequate and we 
obtain some information about angles 5 and k. Figs 5 and 6 show calculated and observed profiles 
for three well known 7-ray pulsars. The energy interval for electrons (or positron) were chosen for 
following reason. Most pulsars exhibit two sharp peaks within a period if the energy of photons lies 
in the region > 1 MeV. Knowing the critical frequency of synchrotron radiation ~ eH^/mc 
we can estimate the Lorentz-factor of particles 7 2 > mcE^/eHh. For the magnetic field of order 
H ~ 10 12 Oe we obtain 703IO or for the particle velocity (3 > 0.99. Specific values of (3 for Figs 7 
were selected in a way to obtain the best fit of peaks width. It is obvious that the greater is the 
particle energy the narrower are the peaks. 

It is believed that the pulsar in the Crab Nebula has an inclination angle close to tt/2. That 
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Fig. 2. — The radiation phase profiles for (3 = 0.9 and angles: a — k = 90°; b — k = 30°; c - 
k = 6° 
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Fig. 4. — The profiles of radiation constructed via exact equation (3) (dashed line) and approxi- 
mation (10) (solid line) for (5 = 0.999 and angles k = 30°, <5 = 81° 
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is why we choose k = 85°. Then the best fit to the observed pulse profile measured by BATSE 
(Fishman) gives the value 5 = 15° see Fig. 5. The next figure shows the Vela pulse profile for 
energy of photons above 100 MeV detected by EGRET (Fierro et al.) along with calculated profile 
for angles 8 = 10°, k = 88° and = 0.999. The difference occurs in a bridge between the two 
pulses. It can be caused by the dispersion of particle orbits in the pulsar magnitosphere. Next 
two figures show Geminga profiles. The first one of Fig 6 displays observed 7-profile for E 1 > 100 
MeV (Daugherty and Harding) along with calculated one for 5 = 89°, k = 10° and (3 = 0.992. Soft 
X-ray profile is plotted in the second of Fig. 6 in the 0.53-1.5 keV band measured by ROSAT 
(Halpern and Ruderman) and a profile calculated for angles 5 = 80°, k = 4° and (3 = 0.94. We see 
good agreement between the model and observed profile though it is thought that the soft X-ray 
radiation is due to thermal black body emission (Halpern and Ruderman). 



The profiles shown on Figs 5 and 6 are calculated for moderately relativistic electrons. In doing 
this we tried to fit the shape of profiles. But all the calculations were made for a single particle. It 
is evident that in a real situation radiation is emitted by an ensemble of particles with a spread in 
energies and orbit orientation. It is evident that on averaging over this parameters the pulse width 
can grow significantly. But without regard on averaging we can obtain some information about 
the pulsar magnetic axis orientation from the relative position of pulses in the radiation profile. 
According to equation (5) the maximum of radiation occurs at given by 



If we know the phase distance between the pulses, we will find the relationship between angles 
S and k by use of equation (11). Here we suppose that 7-radiation of pulsars is a synchrotron 
radiation of particles with their orbits orthogonal to the magnetic axis of pulsar. The probability 
of detection of such pulsar in a 7-ray band according to this model is equal to p = sin/t. Thus, 
with the most probability we can detect pulsars with its magnetic axis orthogonal to the axis of 
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Discussion 



$7ti 5 2 = ± arccos(cot S cot k) 



(11) 
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Fig. 5. — Profiles of 7-ray pulsars Crab and Vela. Smooth lines show the model results, ghis- 
togramms display observed profiles. 
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Geminga Geminga 
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Fig. 6. — Profiles of Geminga pulsars. Smooth lines show the model results, ghistogramms display 
observed profiles, left — 7-ray profile; right — X-ray profile 

rotation. Thereby the gap between pules in 7-ray region should be close to 0.5 of period. That is 
what we see in most cases of detected 7-ray pulsars. 

The discussed model does not pretend to give entire interpretation of the pulsar radiation. 
The main purpose of this paper is to show how to find the light profile of a rotating source with 
an indicatrix of radiation which is given in a rest frame. We chose the synchrotron radiation 
mechanism as an widely spread example. Considered method can be used for curvature radiation 
or more complicated source of radiation. 

This work was supported by Ministry of General and Professional Education of Russia, Com- 
petitive Center of Fundamental Research. 
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